On 
On 
On 

3 



in 



oo 



3 



Ref. SISSA 123/98/FM 

Reduction of bihamiltonian systems and 

separation of variables: an example from 

the Boussinesq hierarchy 

G. Falqui 1 , F. Magri 2 , and G. Tondo 3 



1 SISSA, Via Beirut 2/4, 1-34014 Trieste, Italy 

E-mail: falqui@sissa.it 

3 Dipartimento di Matematica, Universita di Milano 

Via C. Saldini 50, 1-20133 Milano, Italy 

E-mail: magri@vmimat.mat.unimi.it 

Q\ ■ 3 Dipartimento di Scienze Matematiche, Universita di Trieste 

Piazzale Europa 1, 1-34127 Trieste, Italy 
\^q ' E-mail: tondo@univ.trieste.it 

o 
o\ 

-f— > ■ 

a- 

Abstract 

We discuss the Boussinesq system with £5 stationary, within a general framework 
for the analysis of stationary flows of n-Gel'fand Dickey hierarchies. We show how a 
careful use of its bihamiltonian structure can be used to provide a set of separation 
coordinates for the corresponding Hamilton-Jacobi equations. 
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1 Introduction 

It is well known that the Boussinesq equation uu = o(—u xxxx + Au x + 4uu xx ) represented 
as a first order system of equations in the variables (u, v) 

uti = ~u xx + 2v x 

vt 2 = Vxx ~ 2/3 u xxx + 2/3 uu x 

is a member of an infinite hierarchy of evolution equations. These equations are defined as 
follows. In the space M.% of the third-order differential operators L = d x + u{x)d x + v(x), 
one considers the third root La of L, i.e. a pseudodifferential operator of the form 

CO 

L3 = 8 + Y J Qi( u ^)d- i (1.2) 

satisfying (X») 3 = L. By means of the powers of such a third root, one defines an infinite 
family of commutative flows on A4 3 as 

J-L=[L,(L§) + ] (1.3) 

where (■)+ is the projection on the purely differential part of a pseudodifferential operator. 

Stationary flows for such a system correspond to that subspace of operators L satisfying 

the (non-linear) relation [L, L3] = for some p. 

In this paper we will discuss the case p = 5. The corresponding equation of the hierarchy 

reads: 

u ts = -l/9u xxxxx + 5/9 u x u xx + 5/9 uu xxx + 5/3u xx v + 

5/3u x v x — 5/9 u u x — 10/3vv x 
v t5 = -l/9v xxxxx + W/9u xxx v + 20/9u xx v x + 5/3u x v xx + 

5/9 "uw^xa: H — 5/3 Va; — 5/3 ot xx — 10/9nvti x — 5/9 w u^, 



( (5) 



y 



and one can see that the space M3 of the zeroes of — — can parametrized by the space 

of the Cauchy data of two fifth order ODE in two variables. On such a ten-dimensional 
manifold we consider the (restriction of the) first four non-trivial flows of the Boussinesq 
hierarchy 

J-L=[L,(L§) + ], p= 1,2,4,7. (1.5) 

Our aims are to study the bihamiltonian structure of these equations, and to show how 
a careful use of such structures explicitly provides coordinates in which the corresponding 
Hamilton-Jacobi equations are separable. 

(5) 

In the first part we show that on M3 one can construct a linear pencil of Poisson 
brackets {F, G}\ = {F, G}\ — \{F, G}o with the following properties: 

1. they have two Casimir functions H(X) = Hi \ + H 2 and K{\) = H 3 A 3 + ^ X 2 + H 5 A + 
Hq] the functions H\ and H% are Casimir functions of {F, G}o, the functions Hi and 
Hq are Casimir functions of {F, G}i; 



2. the functions H2,H^,H^ and Hq are the Hamiltonian functions, with respect to the 
Poisson bracket {F, G}o, of the four vector fields under scrutiny, while the functions 
Hi, H3, H^, H§ are their Hamiltonian functions w.r.t the Poisson bracket {F, G}\\ 

3. all these functions are in involution w.r.t. the whole Poisson pencil. 



Then we prove that the bihamiltonian structure of the vector fields (1.5) can be used to 
integrate the equations of motion according to the following two-step scheme. 

At first one considers a level surface Sq of the functions H\ and #3. It is an eight- 
dimensional symplectic leaf of the Poisson bracket {F, G}q. The given vector fields are 
tangent to Sq, so they can be restricted to this manifold. Then one remarks that both 
Poisson brackets {F, G}o and {F, G}± reduce to So by a Marsden-Ratiu reduction procedure. 
So one can conclude that So inherits a Poisson-Nijenhuis structure, i.e., So is a manifold 
endowed with a symplectic form u (induced by {F, G}o) and with a compatible Nijenhuis 
tensor N (induced by {F, G}i). 

The second step is to construct a set of separation variables defined on So. To this end, 
we shall use a kind of coordinates introduced in [18] under the name of Darboux-Nijenhuis 



coordinates. They were exploited in [20], as separation variables for the Henon-Heiles type 
systems obtained from the stationary flows of the KdV hierarchy. In this paper, we explicitly 
compute such coordinates for the stationary Boussinesq flows, and finally we show that they 
are separation variables for the Hamilton-Jacobi equations associated with the Hamiltonian 
functions (H2, H^, H^,,Hq). 

In our opinion these flows provide a good example of the tight connection between 
separable coordinates and the bihamiltonian structure of soliton equations. We have chosen 
to study in detail this problem, in order to show that these methods can be effectively used 
to treat systems that are more involved than the KdV cases (see, e.g., |l], 0, ||, 13]), which 



fall into the class of bihamiltonian manifold of maximal rank discussed in [|16| . Accordingly, 
we will illustrate and display the appropriate computations, rather than proving general 
propositions. However, we would like to point out that these kind of results are not specific 
of the Boussinesq equations. Indeed, as it will be clear from our arguments, they have 
analogues holding, mutatis mutandis, for a wide class of stationary flows of n-Gel'fand- 
Dickey hierarchies. Such a formalization will be the subject of a further publication. 

2 The Central System and its reduction 

To obtain Lax pairs and Hamiltonian structures for the £5 stationary Boussinesq system 
(hereinafter the Bsq§ system), we recall the set up for the KP theory discussed in [0, 14 1. 

We consider the space M of sequences of Laurent series {_£f( fe )} fc=1 __ -00 (the currents of 
the theory) having the form 

H (o) = 1; H {k) = z k + ^ Hfz-\ (2.1) 

i>i 

where Hf are scalars. On A4 we define a family of vector fields as follows. We associate 
with a point {-ff' fc '}fc=o,...,oo i n -M. the linear span H+ of the currents H^ k >: 

H + = (1,H^,HW ...), (2.2) 



which is a subspace in the vector space H. of Laurent series in the formal variable z. Then 
we consider the following invariance relations: 

A +ff 0)W cW+J (2.3) 

as the defining equation for the j th vector field of the family. 
Explicitly we have: 

^_ = H { 3+ k) _ H (j) H (k) + £#*#(*-») + j^HfH^- l \ (2.4) 

^ 1=1 i=i 

which will be called in the sequel the Central System. It has the following properties: 

• commutativity: — — — = 0; 

Otj Otk 

• "exactness": -^-H^ = -^-H^. 

Otk Otj 

The connection of the Central System with the usual formulation of the KP theory (see, 



e.g. [27, 12 1 and the references quoted therein) as a system of Lax evolution equations for 
the Sato pseudodifferential operator Q = d + S£i%^~*' as we ^ as the proof of the two 
properties stated above, can be found in @. Here we simply recall that reductions of KP 
to the n Gel'fand-Dickey case are obtained via the constraint 

ffW = z n (2.5) 

a requirement equivalent to the usual constraint (Q n )_ = 0. In particular the Boussinesq 
hierarchy corresponds to the case n = 3. 

It should be noticed that on the subset S n of M. formed by those points satisfying (|2? 
we can read the n equation of CS as the constraint equation 

n 
H (j+n) = z n H (j) _ J2 HiH^ n - l \ (2.6) 

1=1 

which allows to recursively define the Laurent coefficients of H^',k > n, in terms of the 
coefficients of the first n currents. More formally: 

Proposition 2.1 The submanifold S n is the subset ofTi+ given by the equation 

z n {H + ) C Ti+, (2.7) 

i.e., the set of the points where the operator of multiplication by A = z n leaves the space TC + 
invariant. Hence, S n is generated by the first n currents 

H (o) = i, #(!),... ,#(»-!) (2.8) 

over the space of polynomials in A with coefficients in the set 

{Hl 1 ,Hl,...,Hll 1 i },l i eN. (2.9) 

The reduced equation of motion, CS n are simply the restriction of CS to S n . 



In the case n = 3, thanks to Proposition p.l| the degrees of freedom are thus collected into 
the two series 

oo oo 

H {1) = h = z + J2 h i z ~ i andtf (2) =k = z 2 + ^2kiZ-\ (2.10) 

1=1 1=1 

At this level, we have still a dynamical system with an infinite number of degrees of freedom. 

To obtain finite-dimensional dynamical systems representing the restriction of the flows to 

th 
the stationary manifolds, one can remark that the set Z p of zeroes of the p vector field is 

an invariant submanifold for (|2.4f) as well, so that one can restrict the flows on intersections 



(5) 

of invariant manifolds S n V\Z v . In particular, to study Bsqs we shall consider A4 3 = S 3 DZ^. 
The restricted flows are constructed noticing that the coefficients of the Laurent expansion 
of the equations 

give polynomial equations for the coefficients (hi, kj),i,j > 6 which can be recursively solved 

in terms of the first ten variables {hi, ki,i = 1, . . . ,5}. 

Accordingly, upon that substitution, we obtain a system of polynomial vector fields Xj 

(5) 

defined on the phase space Ai 3 whose first members are displayed in Tables 1 and 2. 

Table 1: The Vector field X\ 



hi 


= h-2h 2 , 






hi 


= -h 3 + hi 2 - k 2 , 






h 2 


= -2h 3 -h 1 2 + k 2 , 






k 2 


= -k 3 + h 2 hi - h 4 - hiki, 






h 3 


= —2 h 2 h\ + ^3 — 2 /14, 






k 


= -k 4 - h\k 2 - h 5 - k\h 2 + h\h 3 , 






hi 


= -2h 5 - h 2 2 - 2/ii/i 3 + k 4 , 






k& 


= —h\k 3 + h 2 h\ — h\h 4 + h\ k\- 


- 2 k\h 3 — k\k 2 - 2 k 2 h 2 , 




h 5 


= 3 &5 — 2 k\k 2 — 2 k 2 h 2 — 2 h 2 h 3 - 


-2k l h 3 + 2h 2 hi 2 + 2h x 2 ki - 


- 6/ii/i4, 


h 


= —3k 2 h 3 — k\h 4 — h\k\ — h 2 k 3 - 


- h 2 hi + h 2 2 hi - kik 3 





3 Lax Representation and Poisson Structures 



Another consequence of Proposition ( |2.1| ) is that, (we are sticking to the case n = 3), for 
each j 7^ 3/ there is a non trivial 3x3 matrix V^(A) polynomially depending on A such 



Table 2: The Vector field X 2 



hi 

k 

h 2 
k 
h 3 
k 

hi 

h 5 
h 



-h 3 + hi - k 2 , 
hi + hiki — h 2 h\ — 2 fc 3 , 
-k 3 + h 2 h\ - h 4 - hiki, 
—k\ — h 2 + h§ — h\k 2 — 2 k± + 2k\h 2 , 
-k 4 - h\k 2 - h 5 - k\h 2 + h\h 3 , 

—3 k*> — k\k 2 + k 2 h 2 — h 2 h\ — h\k 3 + 2 hih^ — h\ k\ + 3 k\h 3 — h 2 h 3 
—h\k 3 + h 2 h\ — h\hi + h\ k\ — 2 /ci/i3 — k\k 2 — 2 /C2/12 j 
h\k\ — h 2 h\ + hih§ — 2hik$ + 2k\hi — Akik 3 — k 2 h 3 — k 2 , 
—3k 2 h 3 — kihi — h\k\ — h 2 k 3 — h 2 h^ + h 2 hi — kik 3 , 
h 2 ki - h 2 + h 2 h 5 - h 2 h\k 2 - 2 h 2 k 4 + k\h 2 - k^ + kih 5 - 
h\k\k 2 — 2k\ki — 3k 3 k 2 



that 



^ + Hl " 



1 
h 
k 



V (J) (A) • 



1 
h 

k 



(3.1) 



The commutativity of CS and the exactness property (|2|) implies the zero curvature (or 
Zakharov-Shabat) representation for the restriction of CS to 53. 



■^-V^(X) - ^-V (j) (A) + [V (fc) (A), V {j) (X)] = 



01 



dt k 



(3.2) 



,(5) 



The step to get a Lax problem for Bsqs on -M 3 is now easy. In fact, on the stationary 

d 
manifold of X5 = — — we immediately obtain from ( p.2| ) that the Lax matrix V^> (A) satisfies 

the equations 



_d_ 



y( 5 )(A) = [V (fc) (A),V (5) (A)]. 



(3.3) 



In other words, on the matrix V^ 5 \ all the flows of M.% are Lax evolution equations 
Explicitly, the Lax matrix V' 5 ^ (computed using Eq. ( |2.6| )) has the following form: 



V( 5 )(A)=A 2 





1 
1 

k 3 

2,1 



+ A 



-ko 











1 


h 2 


ill 





h 3 + k 2 


ki - h 2 


-hi 


-h 









+ 



-h 2 hi + hi -hi + h 3 

3,1 







L 



3,2 




h 2 hi - h^ + k 3 



(3.4) 



where 



2,1 


3.2 




hi k 2 + ki - k\ - hi h 3 + h 5 , 
ki h 2 + h 2 2 - h 5 + 2 ki - h 2 , 



Lr,' =hi k 3 



J o 



2 hi hi + 3 k§ + h 2 hi — 2 ki k 2 



(3.5) 



2kih 3 +hi ki - 2k 2 h 2 + h 2 h 3 . 



One can notice that the map {hi,... , A^} i— > V^ 5 ' is invertible, and thus the Lax equa- 
tions ([0]) represent faithfully Bsqs. 

The second members of the Lax pair for the flows X\ and X 2 are given by 



y« 



o l o 

2hi 1 

A + h 2 + ki hi 



y( 2 ) 







A + h 2 + ki 







1 



-h 3 + 2k 2 X-h 2 + 2ki -hi 



(3.6) 



3.1 The Hamiltonian structures 

Besides their Lax representation, the vector fields of Bsqs have another important property: 
they admit as well a bihamiltonian representation, or a Poisson formulation with parameter. 
We have seen that the Lax formulation of the problem comes from the Central System. The 
Poisson formulation can be gotten as follows. 

As it is well known [23, 26, [L7J], on the space q[[X\] of Laurent polynomials with values 
in a Lie algebra g, there is a family of mutually compatible Poisson tensors, Vi associated 
with a family of classical i?-matrices 



R l (X(\)) = (\ l X(\)) + -(\ l X(\))_ . 
Some of them (including the values Z = 0, 1) restrict to the affine manifold 

71-1 

^ := {X e [[A]] |X(A) = X n A + J2 V^.} 

i=0 



(3.7) 



(3.8) 



where A is a fixed element of q. 

A perhaps less known fact, pointed out in p5| , is that Vq and V\ admit further re- 
ductions, leading to generalized Mumford systems. Actually we want to show that the 
Hamiltonian structures of Bsqs , gotten via the reduction of the central system discussed in 
Section |2|, are a Marsden-Ratiu bihamiltonian reduction of the pair Vq and V\ . This choice 
of Poisson pair is by no means accidental. In fact a bihamiltonian equation on q^ 

X = (V 1 - XVo)VF (3.9) 

imply that the polynomial A(A) evolves according to a Lax equation 

We recall that the MR reduction theorem [|1| considers a Poisson manifold (A4,P), a 
submanifold S ^ M and a distribution D C TA4i s such that E = D n TS is a regular 
foliation with a good quotient J\f = S/E. It states that the Poisson tensor P is reducible 
to J\f if the following hold: 

1. The functions on M invariant along D form a Poisson subalgebra of C°°{M). 

2. P(D°) CTS + D, D° being the annihilator of D in T*M. 

The reduced Poisson tensors can be computed according to the following scheme: 

• We choose a covector v£ € T*J\f. 

• We choose a point s E S on the fiber over n, and we lift v\[ to a covector v^ € T*A4, 
that is an extension of tt*v^[ lying in the annihilator D° of the distribution D. 

• Next we construct the vector field (P M ) s v{^ associated with the lifted covector vf 4 
through the Poisson tensor of M. at the point s. The MR reduction theorem insures 
that the vector (P M ) s v^ A is projectable onto J\f. The projection does not depend 
either on the choice of the particular extension v$ or on the point s on the fiber and 
defines {P N ) n v% . 

The bihamiltonian MR reduction theorem Q considers a manifold (M ,Pq,P\) endowed 
with a pair of compatible Poisson structures, i.e. a bihamiltonian manifold. It is a conse- 
quence of the MR theorem stemming from the observation that a bihamiltonian manifold 
admits a kind of double foliation. This suggests to choose S to be a symplectic leaf of Po 
and D = Pi(KerPo)- Then, provided that S is chosen so that the regularity assumption of 
the MR theorem is satisfied, points 1) and 2) above are consequences of the compatibility 
of the Poisson pair (Pq, Pi). Moreover, it states the reducibility of the whole Poisson pencil 
P\ = P\ — A-Poi an d provides A^ with a bihamiltonian structure which can be computed 
following the procedure sketched above. 

In our Bsqs problem we consider the reduction of the Poisson pencil V\ = V\ — XVq, 
associated with the .R-matrices (p. 7]), defined on the space 3x3 traceless matrices of the 
form X(\) = X 2 A + \x\ + xq, with 



4= 10 

1 



Table 3: The reduced tensor Pq 





















-1 














1 























1 


1 
















1 






























-hi 


-2/n 


— fci - h 2 











-hi 


-2Jbi 


— ki - h 2 


-2k 2 


* 















-2h 2 

~k 2 




-k 2 






Making use of the Killing form on sl(3), and of the residue in A, we can parameterize tangent 
and cotangent vectors as matrix-valued (Laurent) polynomials as well: 



X = Xxi + xq; 



A A 2 



The Poisson tensors Vq and Vi have the explicit form [17]: 

.r | | [xi,-] [A,-] wo / [xi,w } + [A,wi] 



X = V {W) o 



Xl 



[A 



o 



U'l 



and 



X = Vi(W)^ 



xo 



xi 



F0, 




[A 




[A wo] 
[w>o,«o] 



(3.10) 



(3.11) 



A long but straightforward computation shows that the Lax matrix V' ' fits into the scheme 

(5) 

herewith outlined, so that via the bihamiltonian MR reduction process the space M^ is 
endowed with a Poisson pencil P\ = P\ — XPq; in the coordinates (h\, hi, h 2 , k 2 , ■ ■ ■ , h$, ^5) 
Po an d Pi have the form displayed in Tables 3 and 4. 

We remark that, in the case of the ^-stationary Boussinesq hierarchy, this procedure 
yields a pair of Poisson structures which coincide with the ones found in [15]. 

The hamiltonian functions for the flows can in principle be computed looking for the 
Casimir function of the Poisson pencil Pi — APo> according the Gel'fand-Zakharevich theo- 
rem pq| . However, since we are dealing with a Lax problem with spectral parameter, they 



are provided as well by the coefficients of the spectral curve associated with the problem, 
i.c: 

= det(/x - V (5) ) = i? - m(#iA + H 2 ) - (A 5 + H 3 X 3 + H A \ 2 + H 5 X + H 6 ) (3.12) 

The explicit expressions of the Hamiltonian functions are: 

Hi =3k§ — 3kik 2 — 3/02/12, H3 = — 'ih\k 2 + 3fc4 — 3fci , 

H 2 = — 3kc,ki — h$k 2 + h^ki — 2h 2 h\h± — ^4^3 + ^3^2 — /13/15 — 

hik 3 ki + h 2 hik 3 + hihak 2 + h\k 2 2 + 3fci 2 /c 2 - ^4^2 + 2k 2 h 2 ki + 
k 3 2 + h 2 - 2h t 2 k i + 2/i 3 fc 4 - 2h 3 k 1 h 2 + 2h 1 h i k 1 + h x 2 h b . (3.13) 

Hi = — 2k\k^ — 2k 2 hi + k^hg, — 3/11/05 — 2h 2 h^ + 'ih\kik 2 + 
4h 2 hik 2 + h 2 + fci - hih 2 h 3 + h 2 k A + h A h 3 - h 2 k z + 
fci/15 + hih^ki — k\h 2 + fc3&2 ^ /i2^i — h\ k\ + 2h\ /14. 
The functions H$ and i^6 are written in Table 5. 



Table 5: The Hamiltonians H^ and Hq 

H^ =h\k\ k 2 + h^h 3 k\ + 3 k\k 2 h 4 + 2h 2 h\k\ + h\h$h 3 — 3/11/02/04 + 3hikik§ + h\h^k 3 

— 2h 2 k\k 4 + k 3 k 2 h 2 — 3 /i2^i h 4 + k 2 h 2 h 4 — 2 h\k^h 3 — 2 h 2 h\k 2 — 3kihi /14 

O O O O Q 

+ 3 h 2 h\k 5 + h 3 k±k 3 + h 2 h 5 + /i 5 /c 4 - /i 5 /ci + /i 3 /c 2 + /12 /ii - /iifc3 

— fe 2 2 /i3 + /ii 3 £:i 2 + /i 2 A;i 3 + /i2 2 ^4 - fci/i 2 3 - 3 /i 4 fc 5 - 4 fci 2 fe 4 + 2 /ii 2 fe 2 2 
+ 2/11/14 + k\h 2 h^ — h 2 h 3 k 3 — h 2 hAji 3 + 2k 4 + 2k\ — /15 — hh 2 h\k\k 2 

H 6 = - k 2 k 2 h 4 - h x k\k 2 - 3 h 2 hi 2 k 2 2 - 3 £^2/13 - k\ 2 h 2 hi 3 + 2 hk 2 2 h 3 + h 2 k 2 k 3 

O OO/I 

— k 2 h^ki + /11/C2 /14 — /15/C1/C4 + 2 /c2 /i2^3 — 2 /C4/11 /C3 + /ci /i2^5 + «i ^2/12 

— fci k 2 k 3 — 2h\k\h 4 + 2k±h 3 k 2 + k\ h\ h 4 — k\ h 3 k 3 + k 4 k 3 k 2 + 3h^kik^ 
+ 2k^h 3 k 3 — 2/14/13/^1 — ?>h 2 h\kik^ — h\k\k 2 h^ + h 2 h\k 2 h^ + 2h\k\k 2 kA, 

— 2 kik 2 h 2 h4 + 3 h 2 h\k\k 2 + h 2 h\h 3 k\ —2h 2 h\h4 : k 3 — h 2 h\ k 2 h 3 + 3h 2 h± k\h^ 

— 2h\k 2 h 3 k 3 J r2h\h 3 k\ki± J r3h\k 2 hi±h 3 — 3h\ k 2 h 3 k\ — k\h 2 k^ J r3h\ k 2 k§ 

— k 2 h 4 h 5 + k 3 k 2 h 5 + h 2 h\k 3 2 + h^h 2 k 3 - h 5 h 3 k 3 + h 2 2 h 3 k 3 - h 2 h 3 2 k 2 

O Q'3Q/I Q O "^ 

— fci/12 /i5 — &1/12 "1 — /iife ^3 + ^1 ^2^1 — 2h\ k 2 hi + &1 /12/C4 — h\h 3 ki 
+ h 3 k 2 k 3 - 2 h 2 k 2 k\ + 3 fci 3 /c 4 + /i 4 2 fe 3 - h 4 k 3 2 + /ci/i 5 2 - 2 fci/^ 2 + A;i 3 /i 2 2 

— k\ h 2 + h 2 h\k 2 k 4 — k\h 2 h 3 k 3 — h\k\h^h 3 + h\h^k\k 3 + k\h 2 hAjn 3 — k\ 



A direct computation shows that 

Proposition 3.1 The coefficients Hj of the characteristic polynomial of the Lax matrix 
V^(A), the vector fields {X\, X 2 ,X^X-j), and the Poisson tensors Pq and P± of Tables 3 
and 4 fill the bihamiltonian sequences of Figure 1. In particular, the hamiltonian functions 
Hj are mutually in involution w.r.t. the whole Poisson pencil. 



4 PN manifolds and the Sklyanin procedure 

In this Section we show how the bihamiltonian structure of the Bsqs hierarchy can be 
exploited in order to get a set of variables which satisfy Sklyanin's separation equations. 
Actually, this kind of variables arise as a set of coordinates which put in "canonical" form 
a Poisson-Nijenhuis (PN) structure obtained by a further Marsden-Ratiu reduction of the 
degenerate Poisson pencil P\. 

4.1 The reduction of the Poisson pencil 

Let M be a bihamiltonian manifold with degenerate Poisson tensors (Pq,Pi) and let Xj 
be the vector fields of a bihamiltonian sequence on M. A possible way to analyze the 
integrability of Xj is to eliminate the Casimir functions of one Poisson tensor, say Po, by 
fixing the values of its Casimirs. Of course, Po can be restricted to any of its symplectic 
leaf So and the Xj, being tangent to So, can be restricted as well to vector fields Xj. These 
are still Hamiltonian on So, with Hamiltonians Hj given by the restrictions of the original 
Hamiltonians to So- However, in general, the bihamiltonian formulation of the restricted 
problem is lost on So- 

Nonetheless, in our situation a specific feature of the pencil allows us to reduce the 
bihamiltonian structure to any generic leaf So of Po- Indeed, the Poisson pencil P\ induces 



on So a Poisson-Nijhenuis (PN) structure |23|, 22||. With respect to such a PN structure, 



the fields Xj admit a formulation, to be introduced in Section a, powerful enough to allow 
us to find a set of separation variables. 

Lemma 4.1 Let us consider the vector fields 

Zi = ^~, Z 2 = — + 2— , (4.1) 

oks 0K4 a/15 

and let T> be the distribution generated by Z\ andZ 2 . Then it holds that: 

£*iZi+W p i " AP o) = Z 1 AW 1 + Z 2 A W 2 ; (4.2) 

where LxP is the Lie derivative of P along X , <f>\ and <p 2 are generic smooth functions and 
the Wi 's are suitable vector fields (which depend on (j>i and (fi 2 ) 

Figure 1: The Bihamiltonian Sequences 

dHi dH 2 

p o/ \Pi 

X{ 



dH 3 dH 4 dH 5 dH 6 

Po/ \ Pi 

X\ X4 X-j 



10 



This condition ensures that the space of functions annihilated by T> is actually a Poisson 
subalgebra of C°°(M^ ). Moreover T> i; 
we can apply the MR theorem to state 



subalgebra of C°°(M^ ). Moreover T> is generically transversal to the image of Pq. Hence 



Proposition 4.2 The pencil P\ is reducible to a pencil P\ on Sq. Since P$ is invertible 
on So, that is, uj = P ~ is a symplectic two-form, we can conclude that So admits the 
structure of a PN manifold, with respect to the pair (Po, N = PiP ~ )■ 

It is convenient to adapt the coordinates on M to the distribution T> as {hi, . . . , h^,g\ = 
2/04 — /15, v 1 = H±,V2 = H^}. In this way, the reduced Poisson tensors can be obtained from 
Po and Pi (written in the new coordinate system) simply deleting the last two rows and 
columns. They are represented by: 



Po 









0-1 

10 

110 

10 

10 








1 





-1 
-1 





-h 

-3fci + /l 2 2/i 2 



-1 








hi 3/ci — /12 
-2/i 2 




(4.3) 



Pi 



10 













-2/ii 

-h 





hi 

-h 

-ki -h 2 

-k 2 -k 2 

h 3 + hi 

hxh 




hi - 2h 2 

-k 2 

-k 2 



-2h 4 + k 3 + 2hxki - 2h 2 hi 

3kj - h 2 ki 

gi — hik 2 — k\ — h\ 





(4.4) 
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The adjoint of the Nijenhuis tensor is 



N* 



-h 


k 2 


k-2 





N l,5 





N 1 ' 7 


3kik 2 - k 2 h 2 





h 


4fci 





-in 2 + h 3 





N 2 > 7 


AT2,8 


-hi 





-h 2 


k 2 


-h 3 - hx 2 





N 3,7 


AT3,8 








2 hi 


h 








h 3 + /li 2 


AT4,8 














kx 


k 2 


A: 2 




















h 





k 2 


1 











2ht 





h 2 











-1 











-hi 


2h 2 -h 



(4.5) 



where 



N 
N 
N 



1,5 



2,8 



4,8 



-2h 2 hi + 2h A -k 3 - 2hih, 
-6h 2 k\ - h 2 2 + 7k\ 2 — gx, 



N 



N 

.3.8 



1.7 



3h 2 2 



h 2 — gx + hxk 2 + kx 
2h 2 kx + gx - kx 2 , 



-2 h 4 + k 3 + 2 hxkx -2h 2 h 



N 



2,7 



Ahxk 



1KX, 



N 



3,7 



-2h 2 hx 



The separation of variables (SoV) is an outcome of the spectral analysis of N* , and 
relies on the introduction of some special set of coordinates. This will be the subject of the 
next subsection. 

4.2 The Darboux-Nijenhuis and the Hankel— Probenius coordinates 

Let (M,Pq,N = PxPq ) be a PN manifold of dimension 2n. In jig] , it was proved that 
if N has n functionally independent eigenvalues (Ai, . . . , A ra ) (otherwise stated, when N is 
maximal) one can introduce a set of coordinates (A := (Ai, . . . , A n ); /x := (/Ui, . . . , /i n )) such 
that Pq and N take the form 



I 
-I 



N 



A 
A 



(4.6) 



where / is the n x n identity matrix and A := diag(\x, ■ ■ ■ , A n ). It was also remarked that 
while the first n coordinates are the eigenvalues of N the remaining ones can in general be 
constructed by quadratures. Since these coordinates are canonical w.r.t. Po an d diagonalize 
the Nijenhuis tensor they will be referred to as Darboux-Nijenhuis (DN) coordinates. 

Another set of coordinates, denoted as (/ := (/i, ... , f n )', c := (ci,... ,c n )) can be 
fruitfully introduced. The idea is the following: (minus) the coefficients Cj of the minimal 
polynomial of the Nijenhuis tensor 



(Det(A-AO)- 



A'' 



X>A- 



(4.7) 



i=X 
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satisfy the following remarkable recursion property: 

N*dci = dc 2 + c\dc\ 
N*dc 2 = dcs + c 2 dc\ 



(4. 



N*dc n = c n dc\. 

In complete analogy, we seek for a set of complementary variables {/i, . . . , /„} such that 

1. their differentials generate an n-dimensional distribution in T*Sn, complementary to 
the one generated by the dc^s; 

2. they satisfy the same recursion relation: 



N*dh 

N*df 2 



df 2 + adfi 
dh + c 2 dfi 



(4.9) 



N*df n = Cndfi 
A characteristic feature of the coordinates (/; c) is that Pq and N* take the matrix form 



where 



H 



Pa 







H 



-Cl 



H T 



1 

-Cn-l 



N* 



F 
F 



(4.10) 



Cl 


Cl 




C, 


1 















1 






(4.11) 



whence the denomination of Hankel-Frobenius (HF) coordinates. 

It is not difficult to show that the relation between the (/; c) coordinates and the DN 
coordinates (A; fj,) is given by the system 



A n_^ A n-* = -Q( A _ A .) 

i=\ 
n 



i=\ 



(4.12) 



I'k 



1, ■ ■ • ,n). 



i=X 



The advantage of using the Hankel-Frobenius coordinate system is that, in the cases at 
hand, they can be computed in a algebraic fashion from the recurrence ( |4.9[ ). 

As a final remark, we note that, owing to the last equations ( |4.8[ ) and ( [4.9[ ), the Hamil- 
tonian fields X = Pndci and Y = Pndf\ are Pfaffian quasi-bihamiltonian vector fields in 



the terminology of P, 20]. Thus these arguments prove the existence of such fields on any 
maximal PN manifold. 
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4.3 Back to Bsq5 

By computing the minimal polynomial of N* of Eq. ( |4.5| ) one gets the c's coordinates 



c\ = k\ + h 2 , c 2 = 2k\ - h\ - gi + 2h\k 2 - h 2 k\, 

c 3 = -3kf + 2fci/i| + 2k 1 g 1 - 2h\k\k 2 - 2h x h 2 k 2 + 2k 2 h A - k 3 k 2 - h 2 k\, (4.13) 

C4 = k 2 h 3 — g\k\ + h 2 k\ — h\k\ + k\ + 2h 2 h\k\k 2 — 2k\h 2 hi + k\k 2 k 3 — h\k\ . 

The /'s coordinates can be found as follows. Taking into account the homogeneity 
properties of the theory with respect to the weights [A] = 3, [hi] = 1 + i, [hi] = 2 + i, we 
see that [Pq] = —8, and [q] = 3i From ( [4.11| ) we get that the weight of /i must be equal 

to —4. Making the ansatz /i = — one finds, using the recursion relation fl4.9| ) that the f s 

are given by: 

, 1 . ~h 2 f , , gi + fe|-2fcg 

11 = -J- ,J2 — —, , h — —til H ; 

fe 2 A: 2 k 2 

kf + fe 2 fcf - k x h\ - g x k x 

h = h ~ h 4 + h x h 2 + — — 

k 2 

The Hankel-Frobenius coordinates associated with the PN structure allow us to recon- 
struct the spectral curve ( 3.12| ), and the Sklyanin separation equations |2g] as follows. In 



analogy with the corresponding arguments known from the KdV case [||, [l3|], we consider 
equations ( [4.12 ) as defining two fundamental polynomials associated with the problem: 



Si = a4 - E c * A(4 ~ l) ; Q 2 = E /* A(4-i) ; ( 4 - 15 ) 



=i j=i 



Then, a direct computation shows that on the zeroes Xj of Q\ and the corresponding values 
Hj = Q 2 (Xj), the equation of the spectral curve (3.12) is satisfied: 



^ = fiji^Xj + H 2 ) + {X) + H 3 X 3 j + H^X) + H^Xj + H 6 ). (4.16) 

When we restrict the problem to the leaf So, we have to consider Hi = v\ and H 3 = v 2 as 



fixed parameters. Hence, we obtain Sklyanin's separation equation 28], i.e. 4 relations in the 
canonical coordinates (A;/x), depending on the 4 commuting Hamiltonians H 2 ,H4,H^,,Hq. 
As it is well known, they imply that, in analogy with the classical Stackel theory of sepa- 
ration of variables for Hamiltonians depending quadratically on the momenta, a complete 
integral W = W\ + W 2 + W 3 + W± of the Hamilton-Jacobi equations 

dW 
Hi(\,—) = ai, i = 2,4,5,6, 

can be gotten by solving the four decoupled first-order ODE 

(dA 2 ) = ^ 2 ( v i A i + a 2) + ( A i+ u 2A J 3 + a4A J 2 + a 5 A i + a 6 ) j = l,---,4, (4.17) 
a problem which can be treated by means of algebro-geometrical techniques. 
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5 Nijenhuis chains and Levi— Civita separability 

Here we introduce the concept of Nijenhuis chains as a natural extension of the Lenard 
sequences and provide and alternative and somehow more classical proof of the fact that the 
Darboux-Nijenhuis coordinates are separation variables for the Hamilton-Jacobi equation 
corresponding to each Hamiltonian vector field Xj of our problem. Remark that, when 
dealing with algebro-geometrical SoV in Section |], the original Hamiltonians ( [3.13 ) of the 



ten-dimensional "Gel'fand-Zakharevich" problem of Proposition |3.1| were not involved in 
the construction of the DN and HF coordinates, and entered the separability proof only 
at the last stage. The reason was that they, not being invariant under the distribution T>, 
cannot define a Lenard recursion relation on So, at least in the classical form enlightened, 



e.g., in [p3|, 11]. Here, we will consider their restriction Hj to the symplectic leaf Sq in 



connection with the following criterion for separability due to Benenti H: 

Proposition 5.1 Let X = P^dK a Hamiltonian vector field on a symplectic manifold 
(M,uj = P ~ ) and (q;p) (q := (qi, . . . ,q n ); P '■= (pi> • • • iPn)) o, set of Darboux coordinates. 
These coordinates are separation variables for the corresponding Hamilton-Jacobi equation 
iff X admits n integrals of motion {.£Q}i<j<„, ; such that 

*(^)*0, (5.1) 

and which are in separable involution w.r.t. the chart (q;p), i.e. satisfy the n relations 

r K K , ._dK t d Kj dK % dK 3 _ 

where no summation over k is performed. 

In order to exploit this criterion, we introduce the notion of Nijenhuis chain which 
generalizes the classical Lenard recurrence satisfied by the (normalized) traces of powers of 
the Nijenhuis tensor I k = ^TrN k , 

N*dl k = dl k+1 => {N^dh = dl k 



and the ones satisfied by the HF coordinates (/; c) as a consequence of d4.8|j4.9|) . 



Definition 5.2 Let (M,P,N) a PN manifold and {Kj}i<j<n n independent functions 
which enjoy the following property w.r.t. N* 

n 

dKj = Y, a jk (N^dKx j = 1, ■ ■ ■ , n , (5.3) 

fe=l 

where ajk are the entries of an invertible matrix-valued function a. This recurrence relation 
(and the corresponding one Xj = PdKj on the vector fields) will be referred to as a Nijenhuis 
chain with generator K\ [X\ = PdK\). 

Nijenhuis chains were implicitly considered in pjl where it was proved that the functions 
Kj admitting the representation (|5.3j ) are in involution w.r.t. the Poisson tensor P. Here, 
we want to state a stronger property of such functions, namely that they are in separable 
involution w.r.t. the DN coordinates. This can be easily proved, by computing {Ki,Kj}\ k 
in the DN coordinates, and recalling that the Xj are eigenvalues of iV*. Therefore in virtue 
of Proposition |Tl] we can state: 
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Proposition 5.3 The Darboux-Nijenhuis coordinates are separation variables for each func- 
tion Kj belonging to a Nijenhuis chain. 

To prove that Bsqs fulfills Benenti's theorem, it is thus sufficient to verify that the restric- 
tions Hj, of the Hamiltonians Hj to a generic symplectic leaf So, belong to a Nijenhuis 
chain. One can directly check that this is the case: indeed, choosing as generator K\ the 
Hamiltonian 



Ho =fc. 



/i 4 /c 3 + h 3 k\ 2 + h 2 2 h 3 



2h\k2 — 3fci &2 + /14 — hikiks + h%h\k^ + h\k2h^ 



2 /11/12/14 — kxk-ihi — 2 /ci/12/13 + 2 kihihz — k\V2 + &2<7i — fefi + hzgi - h\ g\ 



upon the identifications K2 = H^,K 3 = H^^K^ = Hq, the matrix a of the Nijhenhuis 
recurrence is given (in the coordinates (/; c)) by 



a = 



1 













/1C1 - 


-/ a 


A 








/1C2 - 


-h 


-/ici 


./1 





/1C3 - 


-k 


-C2/1 


-/1C1 


/l 



(5.4) 



6 Summary 

In this paper we have discussed, on the ground of the concrete example of the £5 stationary 
Boussinesq system two related items. In the first part of the paper we have sketched a way 
to obtain stationary reductions, and their Hamiltonian structures, of the n Gel'fand-Dickey 
theories in a systematic way. The key points were the study of the Central System (|2.4j ) 
and of its stationary manifolds, and the zero-curvature and Lax representation ( |3.2D one 
naturally gets for the resulting flows. Hamiltonian structures are induced on the stationary 
manifolds by means of a specific bihamiltonian reduction procedure for a pair of Poisson 
structures, defined on the space of Lie algebra-valued polynomials in the spectral parameter 
A. This procedure, whose justification rests on the link between Lax and Poisson formulation 
for Hamiltonian flows, lead us to frame such theories as the Gel'fand-Zakharevich theories 
of a degenerate pencil of Poisson structures, generically (i.e. for n > 2 of non maximal rank 
(see Section ||)). 

In the second part of the paper (Section |] onwards) we addressed the issue of how one 
can solve these integrable dynamical systems by the method of separation of variables. Our 
main theme was to study some of the relations between the classical Hamilton-Jacobi theory 
(and also in its more recent aspects related to the notion of complete algebraic integrability) 
and the Hamiltonian structure we got in the first part. To this end, we had to abandon the 
GZ theory, and, by means of a further (but no more bihamiltonian ) reduction process, make 
contact with the theory of Poisson-Nijhenhuis manifolds. We studied the properties of two 
remarkable classes of coordinate systems associated with the PN structures: the Darboux- 
Nijhenuis coordinates and the Hankel-Frobenius coordinates. By using these coordinates, 
we proved the separability of the Bsqs flows in two ways. In the first we made implicitly 
contact with the "method of the poles of the Baker-Akhiezer function" of the Russian 
School (see, e. g. [13, ^8|). In the second, via the introduction of the notion of Nijenhuis 
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chain, we proved separability by using an equivalent condition to the classical Levi-Civita 
separability criterion |24|| , 

As a final comment, we notice that the peculiarity of the second reduction, which still 
falls in the class of the MR reductions, is given by the fact that two different geometric 
processes are used simultaneously: the restriction for the vector fields and the projection 
for the bihamiltonian structure. Due to this fact, we were able to retain the bihamiltonian 
structure (in the form of a PN structure) but lost the bihamiltonian formulation for the 
vector fields. This is not a new situation in the Hamiltonian theory of integrable system 
and happens, for instance, for: 

• the integrable Henon-Heiles system and its multidimensional generalizations obtained 
by reduction from the stationary flows of the KdV hierarchy pjl, |3(1 : 

• a large class of potentials recovered from the restricted flows of the coupled KdV 
systems, whose most representative member is the Gamier system 0. 

Such classes of dynamical systems live on a bihamiltonian manifold M of maximal rank, and 
their Nijhenuis formulation, namely a special case called quasi-bihamiltonian formulation |], 
^l| ] can be gotten after a similar reduction. The stationary flows of the Boussinesq hierarchy 
(and of the higher n GD theories) are just examples of bihamiltonian structures of non- 
maximal rank, whose interest may lie in the fact that an analysis similar to the one of [JR]] 
for Poisson pencils of maximal rank is not, to the best of our knowledge, at present available 
in the literature. 
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Table 4: The reduced tensor Pi. The entries P ' are the polynomials in the coordinates 
{h\, ... , k$} written below 






-3 








-3/ii 




2/i 2 - fei 


ki - 2h 2 


pi, 10 











3fci 





&3 


- /ii 2 + k 2 


2h 3 - 2hi 2 - k 2 


p2,10 







6/n 


Ski 


3h 2 


2fe 2 


+ 2/l 3 + hi 2 


k 2 + 4/i 3 + 2hi 2 


p3,10 







3ki 


3k 2 


3A; 2 




p4,8 


p4,9 


p4,10 












-3/ii 2 - 3/i 3 




p5,8 


k 3 — 2h 2 hi — 2h,4 


p5,10 











— 3hiki 




p6,8 


p6,9 


p6,10 




* 











p7,8 



p7,9 
p8,9 




p7,10 
p8,10 

p9,10 





pi, 10 . 


= h 3 - 2k 2 - hi 2 , P 2 ' 10 = -hiki + h 2 hi -h 4 + 2k 3 , 


p3,10 _ 


= k 3 + h 4 - h 2 hi + hiki, P 4 ' 8 = k 3 + /i 4 - h 2 hi + hiki, 


p4,9 _ 


2k 3 + 2/i 4 - 2h 2 hi + 2hih, P 4 ' 10 = -2kih 2 + k 2 + h x k 2 - h 5 + h 2 + 2/c 4 , 


p5,8 _ 


-Zhiki + 2h 2 h x + 2/i 4 - k 3 , P 5 ' 10 = -hih 3 - 5hik 2 - 2k x h 2 - 3k 2 + h 5 + k A 


p6,8 _ 


—6k 2 + kih 2 — 2hik 2 + /c 4 — hih 3 + h$, 


p6,9 _ 


-3ki 2 - kih 2 - Ahik 2 + 2A; 4 - 2h x h 3 + 2h 5l 


p6,10 . 


= —3kih 3 — 8kik 2 — Ak 2 h 2 + hik 3 + h 2 h 3 — 2hih 4 + h 2 hi + 3k*, + hi ki, 


p7,8 _ 


h 2 2 + 2hih 3 -k 4 + 2h 5 , P 7 ' 9 = -h 2 2 - 3hik 2 + Ahih 3 - 3k x 2 + 4/c 4 + h 5 , 


p7,10 _ 


= hik 3 + hih 4 — h 2 hi —hi ki — Ak 2 h 2 — 2kik 2 + 2kih 3 , 


p8,9 _ 


-Ak 2 h 2 - 4/ii/i 4 + 2kih 3 - 3kik 2 + 3k 5 + 2hik 3 - 2h 2 h 3 , 


p8,10 . 


= —2k 2 — /11/15 + 2/ii/c 4 — 2k 2 h 3 — 3kih 4 — 2hiki — h 2 k 3 — h 2 h 4 + 2h 2 hi + 3kik 3 , 


p9,10 . 


= —3k 2 h 3 — kih 4 — hiki — h 2 k 3 — h 2 h 4 + h 2 hi — kik 3 
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